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Abstract 

We consider the processor sharing M/M/l-PS queue which also mod- 
els balking. A customer that arrives and sees n others in the sys- 
tem "balks" (i.e., decides not to enter) with probability 1 — 6„. If 
bn is inversely proportional to n + 1, we obtain explicit expressions 
for a tagged customer's sojourn time distribution. We consider both 
the conditional distribution, conditioned on the number of other cus- 
tomers present when the tagged customer arrives, as well as the uncon- 
ditional distribution. We then evaluate the results in various asymp- 
totic limits. These include large time (tail behavior) and/or large n, 
lightly loaded systems where the arrival rate A — > 0, and heavily loaded 
systems where A — > oo. We find that the asymptotic structure for the 
problem with balking is much different from the standard M/M/l- 
PS queue. We also discuss a perturbation method for deriving the 
asymptotics, which should apply to more general balking functions. 
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1 Introduction 



Balking refers to the phenomenon that customers, when forced to wait for 
service, refuse to join the queue. First introduced by Haight f7], balking can 
be specified by a probability distribution associated with the system state. 
Specifically, a customer that finds n customers in the systems upon arrival, 
balks with probability 1 — and joins the queue with probability hn- Haight 
considered several examples of balking functions, including 6„ = l/(n + 1), 
hn = l{n < K} (with !{■} the indicator function), and 6„ = e"*^"". The 
latter case was also studied in detail by Morse |12]. In this paper we shall 
investigate the effect of balking upon systems with processor sharing (PS). 

A model for the round-robin scheduling mechanisms in time-shared com- 
puter systems, processor sharing was first introduced by Kleinrock [8], and 
refers to the service discipline under which every customer gets a fair share of 
the server. It is by now well known that PS is intimately related to random 
order of service (ROS), which refers to the discipline where customers are cho- 
sen for service at random. First studied by Vaulot [15j and Pollaczek [13] , the 
ROS discipline has a long tradition in queueing theory. Pollaczek obtained 
the Laplace transform of the distribution of the steady-state waiting time 
W in the M/M/ 1-ROS queue, by solving a differential-difference equation. 
In fact, the latter was almost identical to the differential-difference equation 
studied by Coffman, Muntz and Trotter [2j for the M/M/l-PS queue. Indeed, 
by comparing these differential-difference equations, it is readily established 
that (see Cohen [3]) 

Prob[V >t]=C ■ Prob[>V > t] (1.1) 

with V the steady-state sojourn time in the M/M/l-PS queue and C a 
constant. A probabilistic argument based on coupling was given in Borst, 
Boxma, Morrison and Niifiez-Queija [T], and the equivalence result was shown 
to extend to other models as well, including finite capacity queues, repairman 
problems and networks. We shall show that for the M/M/1 queue with 
balking the equivalence result also holds. 

The distribution of W does not have a simple representation. Pollaczek 
[T3] was able to invert the Laplace transform and obtained a rather intricate 
but explicit integral representation for the waiting time distribution. The 
integral, along with the method of steepest descent, allowed Pollaczek to 
derive an intriguing asymptotic expression for the tail distribution Prob[VV > 
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t]. This asymptotic expression was rediscovered by Flatto [3]. Morrison 
[TT] considered the heavy-traffic hmit, where the traffic intensity p — 1, 
and derived asymptotic expansions in powers of 1 — p for the sojourn time 
distributions. The tail resuhs of Pollaczek and Flatto were related to the 
heavy-traffic results in Morrison recently in Zhen and Knessl [18j. 

In this paper we consider the M/M/l-PS (or ROS) queue with balk- 
ing. The sojourn time distribution (or waiting time distribution in the ROS 
model) satisfies a differential-difference equation that differs only slightly 
from the one considered by Pollaczek for the systems without balking. How- 
ever, the analysis, and also the system behavior, changes drastically. We 
shall assume that 6„ = l/(n + 1), so that the non-balking probability ex- 
actly matches the share of a server that the customers gets upon arrival. For 
this choice of 6„ the differential- difference equation allows for an exact and 
asymptotic analysis. We obtain the following results: 

(i) An exact spectral representation for the sojourn time density in terms 
of generalized Laguerre polynomials. 

(ii) An expression for the Laplace-Stieltjes transform of the sojourn time 
density. 

(iii) Asymptotic results for tail probabilities when p is fixed; asymptotics for 
the light-traffic case where p — > 0; and asymptotics for the heavy-traffic 
case where p oo. 

The heavy-traffic asymptotics are derived using a singular perturbation ap- 
proach. It is also explained how this approach might be useful for analyzing 
models with more general balking functions. 

1.1 Equivalence relation 

We denote the sojourn time of a non-balking customer that arrives to a PS 
queue with n other customers competing for service by V^, and the waiting 
time of a non-balking customer that arrives to a ROS queue with n other 
customers waiting for service and one additional customer in service by Wn- 
Then we let 6„ and foJi be the non-balking probabilities in the PS system 
and the ROS system, respectively, when there are n customers in the system 
(including the customer in service). 
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Proposition 1. //6q = 1 and b"^ = bn-i, n = 1, 2, . . then 

Vn = (1.2) 

and 

Prob[V >t] = C ■ Prob[y\; > t] (1.3) 

with 

Proof. Borst, Boxma, Morrison and Niifiez-Queija made the observation 
tliat wlienever a service completion occurs in the PS system, each of the 
customers present is equally likely to be the one that departs due to the 
memoryless property of the exponential distribution. In that respect, the 
pool of customers competing for service under PS behaves exactly as the pool 
of customers waiting for service under ROS. Note that the arrival processes 
in both systems can be coupled due to the assumption that fog = 1 and 
If^ = bn-i, n = 1,2, . . .. A similar coupling argument as in |1] then yields 

m- 

Let Mp and JVr denote the number of customers at arrival epochs in the PS 
system and ROS system, respectively. Then, Prob[A^ = n] = 7rop"6o ■ ■ -bn-i 
and 

22n=0P bo---bn 

Similarly, Prob[A/'r = n] = vrgp^^Q ■ ■ ■ b'^_^ and 

P.ob[W>.l^ S"-»^7L-ty'^t"^"^" . (1.6) 

Z^n=0 P ■ ■ ■ 

Upon comparing (11.51) and (11.61) . and using 6q = 1, b"^ = for n = 1,2,..., 
and Vn = W„, the equivalence relation (II. 3p follows. □ 

The M/M/l/K-PS queue can be viewed as a special case of the M/M/1- 
PS queue with balking by choosing bn = ^iin<K — 1 and otherwise. In 
that case the equivalence relation becomes 

1 1 - p^+i 

Prob[V >t] = -- ^ ■ Prob[W > t], (1.7) 
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which was aheady obtained in [T]. For this M/M/ l/i^-PS queue, Knessl 
[9J uses singular perturbation techniques to construct asymptotic approxi- 
mations to the sojourn time distribution. 

We shall consider the M/M/l-PS queue with balking 6„ = in which 

case 

ProblATp = n] = — p, n = 0, 1, . . . (1.8) 

and ^ 

Prob[V>t] = ^ Prob[yy>t]. (1.9) 



2 Problem statement and summary of results 

We consider a processor sharing M/M/1 queue, which also models balking. 
Customers arrive at rate A and the service rate will be denoted by /x. We 
can clearly scale time so as to have /i = 1, and then the traffic intensity is 
p = A//i = A. We let V„ be the sojourn time of a tagged customer that finds 
n others in the system upon arrival. We then define 

V„(t) = Prob[V„ > t]. (2.1) 

With the PS discipline, each customer receives service at rate 1/n when 
there are n customers in service. When a tagged customer arrives we assume 
that he/she will enter the system with probability 6„, and "balk" with the 
remaining probability 1 — 6„. 

The function Vn(t) satisfies the differential- difference equation 

V;(t) = V„„i(t) - (l+phn) V„(t) +p6„V„+i(t) (2.2) 

n + 1 \ J 

with V„(0) = 1. (There is a slight error in jl4j, where this equation was 
previously given.) It is reasonable to define = 1 and have 6„ a decreasing 
function of n. Here we assume that 

6„ = ^— , n = 0,1,2,3,- •• . (2.3) 
n + 1 

Then we define the sojourn time density Pn{t) by Pn{t) = ^^ni^)^ which 
satisfies: 

Pnit) = ^Pn^lit)-(l + ^)pnit) + ^Pn+l{t), t>0 (2.4) 

n + 1 \ n + 1/ n + 1 
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with the initial condition 



PniO) 



n + V 



n > 0. 



(2.5) 



The above can be obtained by integrating (12.21) from t = to t = oo and 
using V„(0) = 1. 

Clearly (12.31) is a very special case of 6„. We can consider also 6„ = 
a/{n + l), since a can be incorporated into the traffic intensity p. However, 
with even a slight change (such as taking 6„ = l/(n + /3) with /3 7^ 1) it seems 
that the problem is no longer amenable to exact solution. We shall discuss 
an asymptotic approach to solving (12. 2p (cf. section 5), which should work 
also for more general 6„. 

We give below various exact and asymptotic expressions for Pn{t). 

Theorem 2.1. The conditional sojourn time density has the following exact 
expression {spectral representation) : 



Pn{t) = ^ Cmil^rn) (pmin, Urn) c''™ * + ^ Cm(z^m) 0m(n, I'm) & 
m=l 



m=l 



where 



and 



-1 + 


1 r 


2m . 






1 r 


-1 + 


2m - 




Cm(z/) 


m 








= n\(^ 






-P 



— p + a/p^ + 4mp 

— p — a/p^ + 4mp 

m-l 



(2.6) 

(2.7) 
(2.8) 
(2.9) 

(2.10) 



Here L^n\z) is the generalized Laguerre polynomial {see p^). 

If we take the Laplace transform of (12.41) and multiply by n + 1, we have 

ppn+m - [{n + 1)(^^ + 1) + p] p^(e) + n p,,_,{e) = -1, (2.11) 

where Pn{0) = Pnit)e~^^dt. Solving the recurrence equation (12. lip , we 
obtain another exact expression for Pn{t), in terms of its Laplace transform. 
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Theorem 2.2. The Laplace- Stieltjes transform of the conditional sojourn 
time density has the following form: 



Pn{e) = MGnY,jHi + MH^ jGi, (2.12) 

1=0 ' l=n+l 

where 

1 

Gn = G„i9) ^ j^"'' - z)'exp ( - {^^dz, (2.14) 

= H^e) ^ r z- - Y^Jexp ( - Y^)dz, (2.15) 

(i+ey 

and 

r = r(e) = — . 
The first two conditional moments of the sojourn time are 

Mn= / tp^{t)dt = —^ + l, (2.16) 
Jo ^ 

Sn = j t'pn(t)cit = y + (2 + -pjn + -p2 + 3p + 2. (2.17) 

Using (I2.12p . we obtain the foUowing asymptotic expansions for p„(t), 
vahd for p > and n and/or t oo. 

Theorem 2.3. For a fixed p > with n, t oo, the conditional sojourn 
time density has the following asymptotic expansions: 

1. n ^ oo, n/t > 1, 

Pnit) = -- — ^ + ^ + 0{n~^). (2.18) 
n n[n — t) 

2. n/t = 1 + = 1 + 0(ri/2), 

p^{t) ~ -erfc( - ) = / e-^'du. (2.19) 

2n V V2/ ny/n J^a/V^ 
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3. Ao <n/t <1 with Aq = ( - p + ^/J^T4p ) /2, 



p„,{t) ~ 1.^"* + ^^"" f-)"e"-*, (2-20) 



n 



where 



r, = r,(-) ^p_(--l). (2.21) 



4. n/t = Ao + A/Vi, A = 0(1), 

p„(t) ^ v^^ + lzy^ e-^Ao"e-*+^°^erfc|-L= L (2.22) 

5. n/it < Ao, 



Expression (12.231) applies also to t — *• cxo with n = 0(1), and gives the ex- 
ponential decay rate of the density Pnit). We note that the right side of (12.231) 
is precisely the m = 1 term in the first sum in (12.61) . i.e., Ci{vi) (j)i{n, ui) e^^*. 
If we start with a fixed large n and increase time t from t = 0, we traverse 
cases 1-5 in Theorem 2.3 in the order given. The leading term in (12.181) is 
Pn{t) ~ 1/n for t < n which corresponds to a uniform distribution. The 
0(n~^) correction term(s) have a singularity as t | n, which indicates that 
the asymptotics become invalid. We also note that if t = 0, (12.181) becomes 
Pn(0) = 1/n — Xjv? + 0{rr^\ which is just the large n expansion of the 
initial condition p„(0) = l/(ri + 1). As tjn increases through one, there is 
a transition region (cf. (I2.19p ) and then for t/n > 1 (but with t/n < I/Aq) 
the density becomes exponentially small, with a rather intricate dependence 
on the space-time ratio, as given in (12.201) . After another transition region 
where t/n ^ I/Aq (cf. (I2.22p ) the density becomes purely exponential in 
t, which corresponds to the dominant singularity in the Laplace transform 
Pn{0), which occurs at 6* = z/i < 0. 

We next consider a small traffic intensity, p — > O"*". We shall consider the 
time scales t = 0(p-^), t = 0(p-^/2) and t = 0(1). 

Theorem 2.4. For p ^ 0^, the conditional sojourn time density has the 
following asymptotic expansions: 
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1. t^c/p-o{p-'), 

(a) n = x/p = 0{p~^) with x > 

Pnit) ~ 1/n. (2.24) 
(h) n^x/p^ 0{p-^) with x^C + ^y/p, 0(1) 

(2.25) 



Pnit) ~ — erfcf 1= 



(c) n — x/p — 0{p ^) with x < C,, 



Pn{t) 



p"H 



■ X 



-3/2 



/27r(C-x) - ^P 

(d) n = XI ^p = 0{p-^'^) with X>C, 



exp|-[x-C + a;log(C/a:)]}. (2.26) 



Pn{t) 



1 



r(i - ^ 

27r V X2 



r,n—t 



(e) X/^ = C>(p-V2) with X = ( + P^'^ Y, 



Pn{t) 



-n/2 -(1-Vp)terfcf^y 



(f) n = X/^ = 0{p-^/^) with X<C, 



i/2 X/2 -C/2 -(l-^)t 



(g) n = 0{l), 



Pnit) ^ ^P '"'"e'-'-e ^'"e 



p^(t)^^p-n/2e-C/2e-(i-VP)t. 



2. t = = 0(p-V2) anc? n = 0(1), 

p^{t)^e-'p-^'^Q^{u), 

where 



(2.27) 

(2.28) 

(2.29) 
(2.30) 

(2.31) 



Q„(u;) = ^(-1)' 



I ™m-n/2-l 



n! m 



m=l 



+E 



2 m! 



n! m 



m=l 



2 m! 



g-m ^(m-l-n) g-c/V^^^. (2.32) 
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3. t, n = 0{l), 

Pnit)=p^^\t)+ppll\t) + Oip'), (2.33) 

where 

row N 1 1 /" r. / 1 e^* , e-* >A 



and 



p-t 



(n + 2)! ^/ + 2 ^/!Vn + 2 n + 2-/ 

^ 1=0 1=1 



n + 1 

Here Br is a vertical Bromwich contour in the 6 -plane with ^{6) > 0. 



Some of the results in Theorem 2.4 for the time scale t = 0{p^^) can be 
derived as limiting cases of Theorem 2.3. However, this is not the case for 
the time ranges t = 0(p~^/^) and t = 0(1). For n, t = 0(1) the leading 
term in fl2.33p corresponds to the tagged customer entering the system and 
no further arrivals entering during his/her sojourn time. 

The result in fl2.3ip is obtained by letting t = cu/^/p and taking p ^ 
in the exact expression (12.61) . If we let — > oo, the m = 1 term in the first 
summation in (12.321) dominates and this verifies the asymptotic matching 
between the scales t = 0{p~^) and t = 0(p~^/^), for n = 0(1). 

Finally, we consider a large traffic intensity, p ^ oo. The structure of 
Pn(t) is different in two cases. 

Theorem 2.5. For p — > oo, the conditional sojourn time density has the 
following asymptotic expansions: 

1. t = Tp = 0{p) and n = Np = 0{p), 

Pnit) = p-i Po(iV, T) + p-2 Pi(iV, T) + 0(p-3), (2.34) 

where 

p m, T) = = , , ^ (2.35) 

°^ ' ^ N-U {N-1){N-U) ^ ' 

and U = U{N, T) is defined implicitly by 



^ - = l-e^-^. (2.36) 



A^- 1 

If N = 1 we obtain the explicit form Pq{1,T) = 
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2. t = T/p = O(p-i) and n = 0(1), 

Pn{t) ~ ^\l - 0" ^0 (2v^v/-e-log(l-0 ) d^, (2.37) 
where Jo{-) is the Bessel function of the first kind. 

We shall compute the correction term Pi{N, T) in (12.341) in section 5 and 
also give some alternate expressions for the leading term Po{N, T), as infinite 
series. The expression in fl2.34p remains valid for n = 0(1) and t = 0{p), 
as well as t = 0(1) and n = 0{p). For N/T ^ 1 we have U ^ T and then 
Po{N,T) ~ 1/iV which is consistent with p„(0) = l/(n + 1) ~ p'^/N. For 
T/N ^ 1, f/ — ^> — 1 and we obtain Po{N,T) ~ e^^ e^"^ , which is consistent 
with Ci 01 e"^^ for p oo and t = 0{p). Note that ui 1/p from (12^1) . 

We remove the conditioning to get the unconditional sojourn time density 
for the PS model as 

oo „ 
n=0 

Then the density p{t) for the ROS model follows from (11. 9p as 

p(t) = (l-e-OPp,(t). 

Note also that the full density, ^^^^(t), for the ROS model is 6(t) +p{t), 
since there is a non-zero probability that W = 0. The exact representation 
for (t) is as follows. 

Theorem 2.6. The unconditional sojourn time density has the exact expres- 
sion 

oo oo 

Ppsit) = $m(^^n^) c'^™ * + ^ C'^(z7^) $^(z7„) 6^™*, (2.39) 

m=l m=l 

where 

Vz/ + 1/ 

We also give the asymptotic results for Ppg{t) and p{t) for the different 
scales of p and t. 
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Theorem 2.7. The unconditional sojourn time density for the PS model 
and waiting time density for the ROS model have the following asymptotic 
expansions: 



1. p fixed with t — » oo 



"^"^ ^ 1-e-P 2v/pT4 

2. p^O 

(a) t = C/p = 0{p-') 

P,sit)-p-'pit)r^^e-'e-</'e-^'-^^'. 
(h)t = uj/^p = 0{p-'l^) 



where 



(c) t = 0(1) 



Qoiuj) = — e cosh(uj / ^/m) . 

^-^ ml 



(71=1 



P.sit) = e-' l + ^it'-2) + 0{p') 



4 
P 
4 



p{t)=pe-' l + ^(t2-4) + 0(p2) 



3. p ^ oo with t = T p = 0{p) 



1 



(2.40) 



(2.41) 



(2.42) 



(2.43) 
(2.44) 

(2.45) 



For fixed p and large t, we removed the condition by using the expansion 
in the region t/n > I/Aq (i-e-, dOSD) in (1^:^ . thus obtaining (^^. 

For a small traffic intensity p, (12.411) on the t = 0{p^^) scale is the limiting 
case of ^rm . as p ^ 0. For the scale t = 0{p~^/^), we used flOTl) in fl238ll . 
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Since p is small, the n = term dominates, which leads to fl2.42l) . When 
t = 0(1), using ([23SD in and the fact that ~ 1 - p leads to 



p,,{t) = (1 -p) p'^\t)+ppT{t)+ppt\t)+0{p') 



which yields (IXlSD . We note that if we let ^ in ([232]), (12:^21) reduces to 
the leading term in (I2.43p . This indicates that the t = 0(1) scale is a special 
case of the t = 0{p~^^^) scale, for small p. 

In the case p — >■ oo with t = 0(p), we used the leading term in (I2.34p in 
(I2.38P and noticed that the infinite sum concentrates near n = p (i.e., = 1), 
which led to fl2.45p . In fact this result is uniform on both the t = 0{p) and 
t = 0{p^^) time scales, for large p. 



3 Derivations of the exact representations 

We first derive the spectral representation (12.61) of the conditional sojourn 
time density. Consider the equation (12.41) and assume that Pnit) has the form 
Pn{t) = e'^^ (f){n). Then (j){n) satisfies the recurrence equation 

(z/+ l){n + l)(t){n) =n(j){n - 1) - p0(n) +p0(n+ 1). (3.1) 

We define the exponential generating function G{z) by 



oo 

Z 



G{z)=Y,-m- (3.2) 



n=0 



Then by (13. ip . G{z) satisfies 

[(z/ +l)z-p\ G\z) + {v + l + p-z) G{z) = 0. (3.3) 
Here we assumed that n(f){n — 1) is finite as ^ 0. Solving (13. 3p . we have 

where C = G{0) = 0(0). Without loss of generality, we let 0(0) = = 1. 

To avoid (f){n) growing like n! as n — > oo, G{z) must be an entire function 
of z, so that — -Rq ~ 1 must be a non-negative integer. The eigenvalues u 
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thus satisfy the quadratic equation = — m, m = 1,2,..., which leads to 
the two sets of eigenvalues given by fl2.7p and fl2.8p . We denote by (prn{n, Vrn) 
and (f>m{n,iym) the eigenf unctions corresponding to the eigenvalues Um and 
Um, respectively. Then for any eigenvalue z/^, using (13.41) we have 

l-'^^z] exp(^-) (3.5) 

m— 1 / 1 \ / , 1 \ m— 1 oo h 



Thus, from (13. 2p we obtain (f)m{n, Vrn) as 

min(?i,m— 1) , , ^ 



mm(n,m-l) ^ i\ / \ / , 1 N2Z-n 

m — 1\ (n\ [Urn + 1) 



We note that as n — > oo, the / = term dominates and (l>min, Um) is asymp- 
totically given by 

n! (z/^ + l)2(— 

0m(^, I^m) - 



(n-m + 1)! (-p)""^ 



(_p)m-l 



n oo. (3.6) 



Alternately, we can use the generating function of the generalized La- 
guerre polynomial (see [TU] ) 

oo 

(l + w)"e-'"^ = J]4"-")(x)w". (3.7) 

n=0 

Comparing (13.51) to (13.71) we see that in our problem, w = —{vm + ^)z/pi 
X = p/ivm + 1)^ and a = m — 1. Thus, by the definition of G in (13.21) . we 
have another representation for (f)m{n, Vm)-, as 

which is (12.101) with v = v^. By a similar calculation, we find that the 
eigenfunctions 4>m{n, Um), which correspond to the eigenvalues Um, also satisfy 
( 12^01) with u = Um- 
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Thus, we can express the conditional sojourn time density as the spectral 
representation in fl2.6l) . with only the two coefficient sequences Cm(^'m) and 
Cmiym) to be determined. 

To determine these coefficients, we ffist obtain an orthogonality relation 
for the eigenfunctions. Since all of the eigenfunctions satisfy (13. ip . we con- 
sider any two eigenfunctions (pmin,^^) and (pm'in,^^') 7^ m'), which 
satisfy 

[{Um + l)(n + 1) - p] (f)m{n, Um) = n0„(n - 1, Urn) + p (t)m{n + 1, Vm) (3.9) 

and 

[{vm' + ^){n+l)-p\ 0m'(n, z/„/) = n0„/(n-l,z/„/)+p0^/(n+l,z/^/). (3.10) 

We multiply (13. 9p by p"' (f)rn'{n,i'm')/n\ and (I3.10p by p"- (t>rn{n,h'rn)/n\, sub- 
tract one equation from the other, and sum over n > 0. This leads to 



n + 1 

n=0 



^ ^ _|_ _[ 

) j— P" 4>m{n, Vrn) <Pm'{n, Z/„/) = 0. 



Since 7^ T^m') we obtain the orthogonality relation 

— P"^ (l)m{n,Vm)(t)m'{n,Vm') (3.11) 

n=0 

By the spectral representation (12.61) and the initial condition (12. 5p . we 
must have 

oo oo ^ 

Crnil'rn) <Pm{n, Urn) + C'm('An) (pmin^Vrn) = —■ 

^-^ ^-^ n + 1 

m=l m=l 

Using (13. lip , we can easily show that, for any eigenvalue Vm or z/^, 
n I.. \- Y.'^=Qp"'^m{n,Vm)ln\ 

LmKVm) - J2-^^(n + l)pn^lin,Um)/nV ^''^'^ 
Using the generating function (13.50 . the numerator in (13.120 is 

oo ri 

y2^4>min,i^m) = Gmip) = (-z/^)"^"' exp ( ■ (3.13) 
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To determine the denominator in fl3.12p . we let G{z,fi) be a solution of 
fl3.ll) with u = fi, whose generating function is given by 



+ 1 \ ^^^^ z 
G{z, /i) = ( 1 z I exp ( ) , R 



Here /i is not necessarily an eigenvalue, and we assume that R < —1 and 
that i? + 1 is not an integer. Thus, (j){n,n), which corresponds to the non- 
eigenvalue n, can be represented as the Cauchy integral 

where the contour C is a small circle in the complex z-plane centered at the 
origin. Since —R — 1 > and R is not an integer, z = p/(/i + 1) is a branch 
point of the integrand in fl3.14p . We also have the binomial expansion 

Y^'' _^ T(R+l + k) f fi + l ^k , 

p V h nR+i)k\ y p ) 

Then from (13.141) . we find that as n — oo, (j){n, /i) is asymptotically given by 
P) ~ , ,x exp 



T{R+l) V p ) "V(/i + l)2 
n\n^ //i + lx*^ / p 



Since (f){n, p) satisfies (13. ip . we have 

[{p + l){n + 1) - p] (l){n,p) = n(j){n - 1, p) + p(j){n + l,p). (3.16) 

We multiply (13.90 by p^ (j){n, p)/n\ and (13.160 by p^ (j)m{'n,h'rn)/n\, subtract 
one equation from the other, and sum over < n < i^, which yields 

^ n + l 

{Urn - /^) j— P" 0m ^m) (p{n, p) 

n=0 

= ^ [MK + 1, i^m) <PiK, p) - 0^(ir, v^) (P{K + 1, p)] . (3.17) 
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We let oo and use fl3.6p and fl3.15l) . which shows that 

+ 1, v^) <P{K, /i) = 0(K™+^-i), 

and 

(j)miK, Urn) ^ + 1, /i) = 0(7^"^+^). 

Thus, the second term inside the bracket in the right-hand side of (13.171) 
dominates the first and we obtain, after dividing both sides hj Um — p and 
expanding for K ^ oo, 

^ n + 1 
^ — ^ n! 

n=0 

(-l)'^ir'"+^(/i+l)^+l(z/„ + l)2— ^-2 / p 



(z/„-/i)r(/?+i)p— 1 v(/i + i 

Next we let ^ in (13.181) . so that R = R{fi) R{i^m) = — By using 
the Laurent expansion of r(-) near a pole 

r(i?+l)~- asi?-*-m, 

^ ^ {m - 1)\ {R + m) 

and then I'Hopital's rule, we find that 

hm (i.^-/i)r(i?+l) =^ \J = / n (3-19) 

i^^i^^l J (m - l)!i?'(z/„) (m- l)!p(l - z/^) 

Thus, by using (I3.18P and (I3.19P and noting that pUm/i^m + 1)^ = we 
let i^' ^ oo and obtain 

^-^ nl 



.„r-(l-.„)exp(^-^). (3.20) 



which determines the denominator in (13.120 . Using (13.130 and (I3.20p in 
(I3.12p . we obtain (12.90 with v = z/^. By the same argument, we find that 
(12.90 is also true for the eigenvalues z/^. This completes the derivation of 
Theorem 2.1. 
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In the remainder of this section, we use a discrete Green's function to 
derive fl2.12p . Consider the recurrence equation fl2.11l) . The discrete Green's 
function Q{6] n, I) satisfies 

pg{9- n + l,l)- [{n + l){l + 9) + p] g{9; n, I) 

+ng{9;n-l,l) = -5{n,l), (n, / > 0) (3.21) 

where S{n, I) = l{n = 1} is the Kronecker delta. To construct the Green's 
function we need two hnearly independent solutions to 

pG{9; n + l,l)- [{n + l){l + 9) + p] G{9; n, I) + nG{9; n - 1, /) = 0, (3.22) 

which is the homogeneous version of fl3.2ip . 
We seek solutions of fl3.22p of the form 

Gn= z''g{z)dz, 
Jv 

where the function g{z) and the path V of integration in the complex z- 
plane are to be determined. Using the above form in fl3.22p and integrating 
by parts yields 

z-g{z)[{l + 9)z-l\ 

- ^ [(1 + ^^) ^ - l\g\z) +p{z-l) g{z)yz = 0. (3.23) 

The first term represents contributions from the endpoints of the contour V. 

If ( 13.230 is to hold for all n the integrand must vanish, so that g{z) must 
satisfy the differential equation 

[(1 +9)z- l\g\z) +p{z-l) g{z) = 0, (3.24) 

and thus 

If the path of integration V is chosen as the segment [0, 1/(1 + 9)] of the 
real axis, then 03.230 is satisfied for n > 1. Thus, we obtain Gn as in fl2.14p . 
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We note that decays as ri — oo, and by scaling z = {1 — y/n)/{l + 9) 
and using the Laplace method, we find that G„ is asymptotically given by 

Gn ^wi^'^i^ ^1 e-''/', n ^ oo. (3.25) 

However, becomes infinite as n ^ —1, and n Gn~i goes to a nonzero limit 
as n 0. Thus Gn is not an acceptable solution to fl3.22p at = 0. 

To construct a second solution to fl3.22p . we consider another path of the 
integration, the real interval [1/(1 + 6'), oo). Then (13.231) is again satisfied. 
Thus, we have another solution of (13.221) . which is given by (I2.15p . Hn 
is finite as n — > — 1, but grows as n ^ oo. By scaling z = nw = 0{n) in the 
integrand of ( 12.15P we find that Hn grows roughly like n! for n large; more 
precisely 

■ I _|_ ()\n+r+l 



Hn ~ n! n' [-^) > n ^ oo. (3.26) 
Thus, the discrete Green's function can be represented by 

g{e;nj) = ^^^^^^^ ifO<n</, ^^.27) 

which has acceptable behavior both at n = and as n — > oo. Here Qo 
depends only upon 9 and /. 

To determine Qo, we let n = / in (13.221) and use the fact that both Gi and 
Hi satisfy (13.221) with n = I. Then we can infer a simple difference equation 
for the discrete Wronskian Gi Hi+i — Gi+i Hi, whose solution we write as 

l\ 

Gi Hi^i — Gi+i Hi = , (3.28) 
P yi 

where Qi = Gi{0) depends upon 9 only. Then using (I3.27P in (I3.22p with 
n = I shows that Qo and Qi are related by Qo = P^~^ Gi/l^-- 

Letting / — > oo in (13.280 and using the asymptotic results in 03.250 and 
(I3.26p . we determine Qi and then obtain Qq as 

Q ^ P r/e 

/!r(r + l)(l + ^) ■ 
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Then, we multiply ( I3.2ip by the solution pi{9) to fl2.1ip and sum over all 
/ > 0. After some manipulation this yields 

oo 
1=0 

which is equivalent to (12.121) . Taking the inverse Laplace transform gives the 
conditional sojourn time density Pn{t) as the contour integral 

Pnit) = 7^ [ Pn{e)e''de, (3.29) 

where Br is a vertical contour in the complex ^-plane, with 3fJ(^) > 0. The 
form in (I3.29P is more useful than the spectral representation for obtaining 
asymptotic results in various limits, such as ra, t simultaneously large. 

To compute the first two moments M.n and 5„, we could expand Pn(6') 
about 6* = 0, but it is easier to derive simple difference equations for the 
moments directly from (12. 4p . By multiplying (12. 4p by t and integrating from 
t = to t = oo we obtain 

pMn+l - (n+ 1 + p)Mn + nMn^l = -{u + 1). 

This equation has the linear solution in (I2.16p . Similarly, for the second 
moment we obtain 

pSn+i -in + l + p)Sn + nSn-i = -2{n + 1) Mn, (3.30) 
whose solution is given by (I2.17p . 



4 Asymptotic results for fixed p and p ^ 

We first assume that the traffic intensity p is fixed. We sketch the main 
points in deriving Theorem 2.3. We first consider n,t oo with n> t and 
use the result in ( 12.12p . To obtain a two term approximation, we need the 
correction terms in the approximations in (I3.25P and (I3.26p . which are given 
by 

Gr 



r(r+l) + -f^r(r+2)--r(r+3))+0(n-2)l (4.I) 
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and 



\ _|_ Q\ n+r+l 



(4.2) 



From (14.11) and ( 14.21) . we note that the first term in (12.121) dominates the 
second, and thus the Laplace transform is asymptotically given by 



Pn{t) 



1 + e 



p9 



1=0 



l\rl 



1=0 

1 (1 + 9) A 



n/ n 



A 

r H 

. n ie\ 



i + e 



pOn 



{l + 9)--y{l-yrdy 



P 



h + h. 



(4.3) 
(4.4) 



(4.5) 



Here 



A = A{e)^'-^^^[2r-e{r + 2)\ 



and we used the Euler-Maclaurin summation formula to approximate the 
sums in (14. 4p by integrals. By scaling 6 = s/n = 0{l/n) and noting that 

^i±^ = p - 1 + ^{2p^ - 9p + 2) + O(n-) 

and r ~ ps/n, the first term, Ii, in (14.51) becomes 



p-1 



n 



-sy 



dy 



p-l\l-e 



n 



Thus, taking the inverse Laplace transforms of Ji and I2 yields 



1 + 



\ — Heavisidefl 

n J n \ n 



n 



(n > t) 



(4.6) 
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and 

n ^ r p-^y 1 

dy ds 



'I — yy-ps/n 



Jo 1 - y 



dy 



n{n — t) 



^ {n> t). (4.7) 



We note that in deriving fl4.7l) , we changed the order of integration and used 
d9 = n-^ds. Then (gS]) and (gZD lead to (KW . 

This analysis suggests that p„(t) is approximately zero in the range t/n > 
1. We shall show that in this sector the density is exponentially small. Before 
doing this, we first investigate the transition region, where t ^ n. 

Thus, we consider n,t ^ oo with n/t = 1 + At'^/^ = i + 0(t"^/2). We 
can still use (14. 3 p but now scale / = y \^ = 0{y/n), and approximate the 
sum by 

1 r°° 

PnW~7 T^T^T^^r^ (l + e)^^y'dy. (4.8) 

^"V ; (i + 0)n+i(^)'-+i ^ J y y k j 

Scaling 9 = zu/ -Jn = 0{l/^/n), and noting that 

and (1 + 6)'^y ~ e'^y, the inverse Laplace transform leads to 

Pn{t) ~ : / e"*/^e^'/2-^^ / e^ydydw 



n 2Tri 



Br 



r J_ f e^'l^e^-^+y+'/^^^dwdy. 



(4.9) 



Note that in this range of {n,t), 



L-v^=(l-^)^ = -A^~-A = 0(l). 



Then by using the identity 

J_ f ^Co^HC,^^^ ^ \ ( _ ^\ 

2m J Br 2v/^ AC J 
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and noting that A = {n — t)/ ^/t, we explicitly evaluate the integral over zu 
in (HH) and obtain i^lM . 

Now we consider n, t ^ oo with t > n. We rewrite fl2.12p as 



oo 



Pn(e) = MGnJ2ji^' + ^ Yl j{HnGi-HiG^). (4.10) 

1=0 ' l=n+l 

The first sum can be calculated exactly by using fl2.15p and making the 
transformation z = + 9) + y, which yields 



J2liHi = e^nr + l)(-^) . (4.11) 



r+l 



The result in fICTD holds for < and 6 > Op = -1 + [-p + ^p'^ + 4p]/2, 
since r(r + 1) has a simple pole at = 9p. We note that 9p = ui, which is 
the first eigenvalue in (12.71) . The second sum in (I4.10p is negligible in view 
of ([325D, fl3:26|l and flCTD . and the fact < on this scale. Using fl213|l . 
(I3.25P and (14. lip in the first sum of (14.100 . then taking the inverse Laplace 
transform, we have 

PnW-TT^ / me'^^'^dO, (4.12) 

2m J Br 

where f{9) = 9 — log(l + 9)n/t and 

r(r + 1) e'^ 



h{9) 



[1 + 9) {-9Y+^n'+^' 



For t ^ oo and n/t fixed we evaluate (I4.12p by the saddle point method. 
There is a saddle point aX9 = 9s = n/t — 1<Q, which satisfies f'{9) = 0. 
Hence, using the saddle point method gives 

„ (t) ^(^.) _ ^tm) 

''''' y^2n tn9s) 
and this leads to (I2.20p . where 

r =r(9) = ^^ = p-(--l) 
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This analysis indicates that (12.201) only holds for n, t — oo with n/t < 1 
and 9s > 9p, so that n/t = l + 9s>l + 9p = Aq. 

There is a transition region where n/t = Aq + A/ \/t with A = 0(1). We 
still use (I4.12P and note that h(9) has a simple pole at 6* = 6*^, and the saddle 
point 9s of f{9) is now close to 9p. We expand the integrand in (14.121) about 
9 = 9p using 

hi9)r.y^i±Iz^e-^ ' 



2^^T4: 9-9p 
and 

/(») = nOp) + fmo- Op) + \nep){e + ■■■ 

- (A„-l)-?.og(A„)-A^ + ^(«-<),)V.... 

Here we also used 1 + = Aq and n/t Aq. By scaling 9 — 9p = S/\/t, 
(I4.12P asymptotically becomes 



Ao 2Ao 



dS, 



where 3?(S') > on the contour Br. Then we use the identity 



2vri J Br S y/27l Ja/^/b 

with A = A/Ao and i? = Aq \ to eventually obtain (E^D- 

Finally, for the scale n, t ^ oo with n/t < Aq, the pole aX 9 = 9p domi- 
nates the asymptotic behavior of Pnif), and (12.231) is obtained by evaluating 
the residue at the dominant pole in (I4.12p . This concludes the derivation of 
Theorem 2.3. 

If we consider a small traffic intensity, p — 0, and scale t = 0{p^^), we 
obtain seven different asymptotic expressions based on different scalings of 
the space variable n, which are given in case 1 of Theorem 2.4. All seven 
results can be obtained as limiting cases of the results in Theorem 2.3, by 
letting p ^ and scaling n appropriately. We omit the derivations here. 
Note that as p the transition line n/t = Aq becomes close to the t-axis. 
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For the time scale t = 0{p ^/^), from the spectral representation we note 
that as p the eigenvalues are 

i/^ = -1 + ^ + 0(p), V^ = -l-^ + 0{p). 
\/m Jm 



Then the eigenfunctions are asymptotically given by 

and 4>m{n,Um) ~ (— 1)*^ 0m(^, '^m)- Thus, all the eigenvalues contribute to 
Pnit) on the scale t = 0{p^^^'^) and n = 0(1), and we obtain (12.311) . 

Now we consider the scale n, t = 0(1) with p ^ and use the result in 
fl2:T2|) . Since r = 0(p), G„ in (imj) becomes 

fiTe 1 

G„~ / z''dz = - — — . (4.13) 

Jo + 1) (1 + ^^)n+l ^ ^ 

By scaling z = jp^gj^ in (12.151) . Hn is asymptotically given by 

1 

pn+l (^f^ 

- n'.^'^T\ (4.14) 

We also have M ~ p/(l + 0). Using (H:T3il and (Hlil) in (l232ll . we find that 
the first sum dominates the second and we obtain 



Hr, ~ r-r^^ / s" exp ( - — — ] ds 



n 



P 

1=0 

n 



(n + 1) (1 + 0)"+2 



1 1 



n + 1 9 



1- ' 



1 + ^)"+i_ 

Then we invert the Laplace transform over time, letting 6 = w — 1, which 
gives 

e"* 1 r 1- w""-^ 



Pnit) ~ — — I e'"'dw 

n + 1 27ri J w — 1 
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n + _ 



This is the leading term Pn\t) in (12.331) . 

To obtain the second term pn\t), we need the correction terms of the 
asymptotic expansions in (14.131) and (I4.14p . It is much easier, however, to 
use a perturbation method to obtain pn\t). Assume that the conditional 
sojourn time density has an expansion in powers of p, as in (12.331) . By using 
the recurrence equation (12.41) . the leading term Pn\t) satisfies 

with the initial condition Pn\t) = l/(n + 1). This can be easily solved to 
regain (14.151) . The second term Pn\t) satisfies 

^P"^^^) - Ji) (t) - „(i)ft) + ^_ „W (t) - ^— nW(t) 

with the initial condition pn\o) = 0. We set 

Pn\t) = ^Vn{t) (4.16) 

and take the Laplace transform of Vn{t) over time, with P„(s) = Vn{t) C^dt. 

After some simplification, we find that Vn{s) satisfies the following difference 
equation: 

( \ -B f\ (n + l)s-"-^-(n + 2)g-"-^ + l 

sVn{s) -Vn-lis) = , n>l 

[n + l)[n + 2)[s - 1) 

with 

^ / N 1 — S 

After some calculation, we obtain P„(s) as 

^ . ^ 1 (n+ l)s-" 



s-1 (n + 2)s(s-l) (n + 2)s2(s-l) 



+ 



^ n ^ n— 1 

n+3 / I 9 ~ 



J-n-1 



/ + 2 ^ 1 + 2 

1=0 1=0 



Then by taking the inverse Laplace transform, and using the relation (I4.16p . 
we obtain Pn\t) in (I2.33p . 
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5 Asymptotic results for p ^ oo 

We shall use a singular perturbation approach to derive the asymptotic ap- 
proximations for large traffic intensities, p — > oo. This method should be 
useful for analyzing models with more general balking probabilities. We 
shall also sketch how the asymptotic results in Theorem 2.5 follow from the 
exact representations in Theorems 2.1 and 2.2. 

We ffist consider the scale t = Tp = 0{p) and n = Np = 0{p), and 
expand Pn{t) in powers of p~^, as in fl2.34p . Using fl2.34p in the recurrence 
equation (12. 4p . the leading term Pq{N,T) satisfies 

dPo N-ldPo_ 

dT N dN N ° ^ ■ ^ 

with the initial condition Po{N,T) = 1/N. We solve this ffist order PDE by 
the method of characteristics. The family of characteristics is given by 

T = + log |1 - A^l + constant, 

where the constant indexes the family. The characteristic T = N +\og{l — N) 
goes through the origin (A^, T) = (0, 0), along the parabola T = —N'^/2. The 
general solution to (15. ip is 

Po(iV, T) = ^^((iV - 1) e^-^). (5.2) 

Using the initial condition in (15. 2p . we determine the function JF(-) from 

N-1 



^((Ar-l)e^) 



N 

If we denote by A^* = A^* ( A^, T) the solution to 

(Ar,-l)e^* = (Ar-l)e^-^, (5.3) 

Pq in (15. 2p becomes 

1 A^ - 1 

PJN, T) = —— — -. (5.4) 

' ' N-1 N, ^ ' 

Setting = N — U in (15. 3p leads to (I2.36p . Thus, (15.40 can be rewritten as 

dOSD. 



27 



Alternately, we can rewrite fl2.35p more explicitly, in terms of an infinite 
series. From (^M), we let U = N - 1 + Uo, where Uo = Uo{N,T) = 
[l — N) e^~^. Then Uq can be expressed in terms of the Lambert W-function 
(see [4]), which satisfies 

e-^o Uo = {l- N) e^"^-i = z. 

We use the series expansion of the Lambert W-function to obtain Uq as 



ml 

m=l 



where the series converges for \z\ < e ^. Thus, U has the following series 
expansion 



U{N, T) = AT - 1 + V — — (1 - AT)™ e"^(^-^-^), (5.5) 



771=1 



which converges for |1 — A^| < 1. The series is always convergent for 

A^ < 1, but diverges for A^ > 1, if T < A^ + log(A^- 1). For example, if T = 
the series converges only for N < Nc = 1.2784, where (Nc — l) e^" = 1. Using 
(15.51) in (12.351) . we have an alternate series expression for Po{N,T): 



oo 



J2 (1 - ^)™"' e'"^^-^-!) /m! 
Po(Ar, T) = . (5.6) 



m=l 



Now we sketch how to compute the correction term Pi{N, T), which sat- 
isfies the following PDE 

dl\ ^ 1-N dPi_]_ N + 1 d^Po N-l dPo 1 , . 

dT N dN N 2N dN^ m dN ^ ' ^ 

with the initial condition Pi(A^, 0) = —l/N"^. This follows from expanding 
Pn(0) = l/(ra + 1) = l/(A^p + 1) in powers of We make the substitution 

Pi{N,T) = -1-P,(N,T) + P*{N,T), (5.8) 
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so that P*{N, 0) = 0. 

We introduce an operator V, which is defined by 



^ d N -1 d 1 , , 

Then by (15. ip we have VPq = and 

Applying the operator V to (15 .Sp . and using (15.70 and (15.100 . we obtain 

We change variables in (15.110 from (A^, T) to (^, r^), where N = ^ + 1] and 
U = Tj, with [/ = U{N,T) given by (I2.36p . Thus, by the chain rule, we have 



dN V 9A^/ 9^ SAT dri 
and 

By implicitly differentiating (I2.36p with respect to A^ and T we obtain 

dU _ T) dU 

ON ~ ^i^ + r] -1)' df~ e ■ 

Then by (15. 9p . after changing variables, the operator V can be rewritten as 

P=i±l^# + ^. (5.12) 

^ + r] drj ^ + 7] 

If we denote Pi{N, T) by -F(^, rj), then by the chain rule we have 

dPo ( dU\dF dUdF ^ , , 

l-^Kr + T777^ = ^i(^'^) (5-13) 



ON V aA^/ 5C aAT d-q 

and 



a^Po A 9f/x dF, dU OF, _ p . ^ 
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Thus, by using fl5.11l) - fl5.14l) and noting that N = + r], we have 



VF 



^ + 1] drj ^ + 7] 



1 



2{^+v) i^+vy i^+vr ■ 

Multiplying the above hj C, + 1]{= N) and using fl5.12p we obtain 



drj 



(e + r/-l)F 



2 ^ ^ + rj ' 



F. (5.15) 



Solving (15.151) with the help of the symbolic computation program MAPLE, 
then imposing the initial condition Pj'(A^, 0) = F(^, 0) = and using (15. Sp . 
we obtain the correction term Pi{N,T) in (I2.34p as follows 



PiiN,T) 



(^-l)(2e-3) 3(e-l)(2e^-2e-3) 
2^ 2^(^ + 7^-1) 

(e - iyi2e + 2e - 5^ - 15) (^ - 1)^(2^^ +4^ + 3) 



2^(^ + 77-1)2 



2e(e + r/-l)= 



2(e-l)(2e-3) 



log 



e + r/-l 



e(e+r/-i) 

where = N — U, rj = U. 

Now we consider some special cases. If = 1, then U 
Thus, ,^ ~ 1, 77 ~ and 

-T 



(5.16) 



by fl236|) . 



Then (15.160 reduces to the explicit result 

9^ 



Pi(l,T)= (2T--)e-^ + 8, 



9 

- e 
2 



We already showed that Po{l,T) = . We relate the explicit result along 
= 1 (n = p) to the spectral expansion (12.61) . From (12.61) we can easily 
show that only the eigenvalues i^i, 1^2 and z/3 are 0{p~^) or 0(p~^) in this 
limit (the others are o{p~'^)). Expanding (12. 60 - 02. lOp for p — 00 and A^ = 1, 
we obtain 



Ci{yi) 01 (p, vi) e 



Ul t 



"1 „~T 



p ^ +p" 
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and 




3T 



Thus (12. 6p agrees precisely with p ^ Pq{1,T) + p ^ Pi(l, T). 

If = and T 00, then f/ —1, 77 ~ —1, and ^ ~ 1 — e^^ e^"^. Hence, 
from dESD, Po(0, T) ~ e"^ e"^ and Pi(0, T) ~ (2r - 3) 6"^ e"^. From (1231) 
we obtain for n = 0, t = Tp —>■ 00 

Ci(z/i) 0i(O, z/i) e^^ * ~ p-^ e-^ e'^ + p-^ (2T - 3) e'^ . 

Again, this agrees with (I5.16P and shows that for = and T ^ 1 (t ^ p) 
only the first eigenvalue ui contributes to the expansion of Pn{t)- 

Next, we consider short time scales, with n = 0(1) and t = r/ p = 
0{p~^). Expanding the conditional sojourn time density p„(t) in the form 
Pn(t) = Qn{T) + 0{p~^) and using equation ( 12. 4p . we have 



with the initial condition Qn{0) = l/{n + 1). Taking the Laplace trans- 
form over the time variable r with Qn{s) = Qn,(T) e""^^ dr, we obtain the 
following difference equation for Qn{s): 



Q'nir) 



1 



Qn+i(r) - Q„(r) , 



n + 1 L 



Qn+iis) - [{n + 1) s + 1] Q„(s) = -1. 



(5.17) 



Solving (I5.17P yields 




(5.18) 



By the inverse Laplace transform, we have 




(5.19) 



Using the identity 




X, y >0, 
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we can rewrite fl5.19l) as 



Br ^ Jo 



2m 
1 



1 [ e'/' (1 - 2)^/^ 



{l-zY / e"* ^ — dsdz. 

2m J Br s 

Then by the inverse Laplace transform 

C-'{le^/^)=M2V^\), 
where W = z + log(l — z) < and Jo(-) is the Bessel function, we obtain 

Qn{r) = ^ V - Jo (2v^V-z-log(l-z) ) rf^, 

which leads to (12.371) . By expanding Qnir) for n ^ oo and r — oo, with 
T = 0{n^) we obtain 



poo 

Jo 



+ 2r pVA^2 + 2T' 



(5.20) 



Then we can easily show that p^^Po{N,T), when expanded for {N,T) — >■ 
(0, 0), gives the same result as in (I5.20p . which verifies the matching between 
the long time (T-scale) and short time (r-scale) results. 

We show how to obtain the results in Theorem 2.5 from the exact repre- 
sentations. First consider n, t = 0{p). Then from (12.71) and (12.81) we obtain 
Urn ~ —m/p and = 0{p). Thus on the large time scales t = Tp = 0{p) 
we have e^™* ~ e~"^"^, while e'^™* becomes exponentially small. Thus all of 
the terms in the first sum in (12. 6p contribute to the leading term for Pn{t)- 
We furthermore scale n = Np and use 



p m\ 

and 



n! ( ^ ) ~ (-p)" (1 - iV)'"-i e^'"^. 
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Thus from (12.61) . on the {N,T) scale, we obtain 

Pn{t) P^' y e"^(^-i) (1 - iV)™-! ^ e-'"^. 

ml 

m=l 

We can show that this is equal to p^^Po{N,T) in (15.61) . 

Next we consider n = 0(1) and t = r/ p = 0{p~^). We use the exact 
representation for Pn(t) in Theorem 2.2 and scale 9 = ps — 1 = 0{p) (s > 0). 
Then r ~ 1/s and from (I2.13P we have 

M --. (5.21) 

sr(l + l/s) ^ ' 

By scaling z = w/{ps) = 0{p^^) in (12.141) . G„ becomes 



1 



{psr+w/s V ^ 

n\T{l + 1/s) 



(5.22) 



r(n + 2 + i/s) (ps)"+i+i/*' 

We also find, from (I2.15p . that Hn is asymptotically given by 

POO 

Hnr^ z"+^/'^e~"/^rfz = s"+^+^/^r(n + l + l/s). (5.23) 
Jo 

Now we rewrite (12.121) as 



I °o J 



Pr,ie) = MGnJ2^m + MHnJ2^Gi (5.24) 



/=0 l=n 



and use (I5.21l) - (l5.23p in (I5.24p . We find that the first term in the right-hand 
side of (15.240 is 0{p~^) and the second term is 0{p~^). Thus, the second 
term dominates the first and we have 



oo J 
l=n 



ly^ z_i r(n + l + l/.) 
p^ Til + 2 + l/s) 

1>^^_,_, r(n + l + l/5) 
Pj^o nn + J + 2 + l/sy 
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which corresponds to p in fl5.18l) . Inverting the transform using dO = 

pds leads to the same result we obtained by the perturbation method. 



6 Discussion 



To summarize, we have obtained both exact and asymptotic results for the 
M/M/l-PS model with non-balking probability 6„ = l/(n + 1). We compare 
our results to the standard model, where bn = 1. First, the spectral repre- 
sentation of Pn{t) for the two models is very different as the standard model 
has a purely continuous spectrum (see also Guillemin and Boyer ^) while 
the balking model has a purely discrete one. 

We recently studied (see [T7j) Pn(t) for the standard model asymptotically, 
and found that ifp = A//i<l and ra, t — > oo the asymptotic expansion is 
different according as n/t > 1 — p, n/t ^ 1 — p, < n/t < 1 — p, n = 0(t^/^), 
and n = 0(1). The scale n = 0{t^/^) is important in obtaining the tail of the 
unconditional density, which for the standard PS model has the form (see 
[HI and [5]) 

where = {l — y/pY and ai and 0:2 are constants. In contrast, for the model 
with balking Theorem 2.3 shows that the structure of p„(t) is different in three 
main sectors of the {n,t) plane {n/t > 1, Aq < n/t < 1 and < n/t < Aq), 
with two transition regions connecting them. For t —>■ 00 with < n/t < Aq 
the asymptotics of Pnit) are governed by the eigenvalue with the largest real 
part and we obtain the purely exponential behavior in (12.231) . which leads to 
f l2.40p for the unconditional density (t) . Thus for the model with balking 
the scale n = 0{t^/'^) is absent. 

If p > 1 the standard PS model has an algebraic tail, with Pps{t) ~ 
^-p/(p-i)^ so that the mean sojourn time is finite for p < 2, the second 
moment is finite for p < 3/2, etc. Then the approximation 



Pn{t) 

n 



1+ p-1 - 

n 



p-1 



applies for n and/or t — > 00. This situation is similar to the model with 
balking in the limit p ^ 00. Here the tail will be purely exponential, but for 
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n and/or t ^ oo we have the approximation in fl2.34l) . which is quite unhke 
the three sectors in Theorem 2.3. 

We derived Theorem 2.5 by both a perturbation method and by using the 
exact representations. The former method should also be useful for general 
non-balking functions 6„, provided that we can write p6„ in the form pbn = 
B{en), where e is a small parameter. Thus pbn is a "slowly varying" function 
of n. For example, this would apply to 6„ = e~^" (used by Morse [12]) if 
c is small. This limit would also apply to repairman problems (or finite 
populations queues) where hn = M — n and M is the customer population. 
Then phn = pM{l — n/M) and we would assume that M — oo (thus e = 
M~^) and p — 0, with pM = 0(1). It is likely that the asymptotic structure 
of all of these models is quite different, and the perturbation method should 
clearly show these differences. 
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